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We study the lensing properties of an asymmetric mass distribution and vector field in General- 
ized Einstein-Aether (GEA) theory. As vector field fluctuations are responsible in GEA for seeding 
baryonic structure formation, vector field concentrations can exist independently of baryonic mat- 
ter. Such concentrations would not be expected to be tied to baryonic matter except gravitationally, 
and so, like dark matter halos, would become separated from baryonic matter in interacting systems 
such as the Bullet Cluster. These vector field concentrations cause metric deviations that affect 
weak lensing. Therefore, the distribution of weak lensing deviates from that which would be in- 
ferred from the luminous mass distribution, in a way that numerical calculations demonstrate can 
be consistent with observations. This suggests that MOND-inspired theories can reproduce weak 
lensing observations, but makes clear the price: the existence of a coherent large-scale fluctuation 
of a field(s) weakly tied to the baryonic matter, not completely dissimilar to a dark matter halo. 



PACS numbers: 04.50.Kd 
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INTRODUCTION 

It is well known that General Relativity can not explain 
the full dynamics of our universe and the structures in it 
sourced solely by visible mass. Two kinds of exotic energy 
density, dark matter and dark energy, must be introduced 
to explain respectively the dynamics of structures and the 
accelerating expansion of the universe. 

An alternative to adding new sources of gravity is to 
alter the response of gravity to the known matter sources. 
Modified theories of gravity replacing dark matter as an 
explanation of flat galaxy rotation curves were first pro- 
posed by Milgrom [fl, 0] and by Milgrom and Bekenstein 

@] , and came to be known as Modified Newtonian Dy- 
namics a.k.a. MOND (for a review of MOND see 0). In 
the contemporary incarnation of this phenomenological 
theory, Poisson's equation for the Newtonian potential 
(j), is altered by the inclusion of a scalar function fi: 



V-(M(rW>|)V£) =4*G N p. 



(1) 



The scalar function is chosen to enhance gravitational 
acceleration at large distance, and thus replace galactic 
dark matter. It turns out that galactic dynamics be ex- 
plained using very simple choices of /x, such as 



m(IW|) 



|V<^|/a 



(2) 



Introducing only one new parameter, a universal acceler- 
ation, do, one can fit almost all galactic rotation curves 
assuming very reasonable mass-to-light ratios consistent 
with prior expectations. 

Despite its successes, most scientists did not regard 
MOND as a viable theory. In part this was due to 
difficulties explaining clusters and possibly some small 



scale structures, but more particularly it was due to 
the absence of a consistent covariant action formula- 
tion of the theory. More recently, these theoretical hur- 
dles were overcome by Bekenstein who proposed 
Tensor- Vector-Scalar gravity (TeVeS). Zlosnik et al.^ 
then showed that TeVeS is part of a larger class that 
they call Generalized Einstein- Aether (GEA) theories. In 
these theories, a new field or fields modify the response 
of the metric to the presence of matter. Most generic is 
the existence of a vector field with a non-zero time-like 
background value. Hence the name Einstein-Aether. 

It has been claimed that the discovery of the merging 
galaxy cluster 1E0657-56 the bullet cluster, rules 
out these modified gravity theories and proves the exis- 
tence of dark matter [2]]]. X-ray images of this cluster 
(left hand side of figure [T]) reveal a bullet-like subcluster 
just exiting the core of the main cluster. Weak lensing 
mass contours (right hand side of Figure [T]) show that 
the lensing centers are not at the luminous centers of 
these bright clusters, rather they are (approximately) on 
the line through these X-ray centers but outside them 
(cf. Figure [1]). The dark matter interpretation is that 
the lensing centers (which are less luminous than the 
X-ray centers) are caused by the presence of collision- 
less cold dark matter (CDM). In the aftermath of the 
cluster-cluster collision, the CDM halos have presumably 
continued along un-impeded, while the hot baryonic gas 
has been slowed. 

It is true that MOND is unlikely to explain the sepa- 
ration of the gravity centers from the mass centers. But 
MOND was never to be regarded as a complete the- 
ory, and certainly not as a theory of time-dependent 
mass distributions. Previous authors 2H, 2fi| have taken 
the MOND-derived Newtonian potential, inserted it into 
TeVeS to show that TeVeS too does not explain the ob- 
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FIG. 1: The bullet cluster [2l|: Shown above in the left panel is a color image from the Magellan images of the merging cluster 
1E0657— 558, with the white bar indicating 200 kpc at the distance of the cluster. In the right panel is a 500 ks Chandra image 
of the cluster. Shown in green contours in both panels are the weak lensing k reconstruction with the outer contour level at 
ft = 0.16 and increasing in steps of 0.07. The white contours show the errors on the positions of the k peaks and correspond 
to 68.3%, 95.5%, and 99.7% confidence levels. 



served weak lensing map without an additional dark mass 
component. This has led the Bullet cluster to be widely 
regarded as direct evidence for dark matter and against 
modified gravity. 

But are relativistic versions of MOND truly ruled out 
by the bullet cluster? The new dynamical degrees of 
freedom of these theories arc known to solve other clear 
problems of vanilla-MOND. For example, it was shown 
[lfj that perturbations glows too slowly in a MOND 
universe for the observed large scale structures to have 
emerged from the fluctuations present at the time of re- 
combination. However, it has been shown in both TeVeS 



111 [12| and GEA growing modes of the new dynamical 
fields introduced in these modified gravity theories suc- 
cessfully seed the growth of baryonic structures after the 
end of silk damping. This implies that observed bary- 
onic structures - galaxies, clusters, etc. - evolved from 
non-baryonic seeds which were not sourced by baryons. 
If these seeds persist, then they might be expected to 
behave in very much the same way as a dark matter 
halo when the baryonic structures with which they are 
associated collide - travelling on, in a straight line unim- 
peded. These non-baryonic structures may then be the 
non-luminous weak lensing centers observed in the bullet 
cluster. Of course, this would mean that even in modi- 
fied gravity theories baryonic structures are surrounded 
by non-baryonic halos - but halos of weakly interacting 
classical fields rather than halos of weakly interacting 
non-relativistic particles. 

MOG is a possible modified gravity theory to 
explain bullet cluster 13]. For our calculations, we 
will focus on the Generalized Einstein- Aether theory. It 
is a general form of "Einstern- Aether" theory 
|l4L Il5l |. In this theory, general relativity is modified 
only by the addition of a dynamical vector field A a , with 
a time-like classical value A a = (1,0,0,0). We mod- 



ify the conventional ansatz for A a in the presence of a 
matter source to allow for first-order perturbations in 
the spacelike, rather than just the time-like components. 
The spacelike component is not sourced by the matter, 
but could emerge from the intrinsic vector growing mode. 
The vector field enhances weak lensing by different fac- 
tors at different locations, with a resulting possible disso- 
ciation of the luminous and gravitational-lensing centers. 

We organize the paper as follows. We first review 
briefly the Generalized Einstein- Aether theory and cal- 
culate the solution for the vector field to first order. We 
then discuss appropriate boundary conditions for the bul- 
let cluster, and perform suitable numerical calculations 
of the resulting metric and convergence. 



THEORY 

The action of the Generalized Einstein- Aether theory 
can be written in the form [l[ 



S 



d 4 x. 



R 



16-kG 



N 



£(9, A) 



Sm- 



(3) 



Here, g is the usual metric (with signature (-,+,+,+)), R 
the Ricci scalar of that metric, and Sm the matter action. 
A is a new dynamical vector field. As shown in [H, @], 
we can obtain a MOND-type limit in the appropriate 
low-acceleration regime by choosing 



M 2 



16ttG n 
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\{A a A 



(4) 



where 



K. = M- 2 K ap ^ a A~<VpA< 7 and 



K. 



a/3 



75 



c 3 6^ 



(5) 
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The Ci are dimensionless constants, whereas M has the 
dimension of mass. A is a non-dynamical Lagrange- 
multiplier field with dimensions of mass-squared, that 
enforces that A a is unit-timelike. The gravitational field 
equations obtained by varying g a @ @, [jj) are 



G a p — T a p + 87rGArT a/3 



(6) 



where the stress-energy tensor for the vector field is given 
by 



1 



-T'Y< 



(a/3) 



-g a pM 2 T + XA a A 



(7) 



where 



T = die and 

(/C Q/3 ^ 



(8) 



Brackets around indices denote symmetrization. Y a p is 
the functional derivative 



Y 



a/3 



= V.A^A^^ 



7 ^ 6g a P 
= -ciKVMiVA^-iVMiVpA")} 

The equations obtained by varying A@ are 

V a {T l J a l3 ) + T' V p = 2XAp, 

with 



(9) 



(10) 



SAP 



For the choice of K. given by equation [H yp = 0. 
Variations of A fix 



A^A, 



-1. 



(11) 



(12) 



The MOND type solutions of the Generalized Einstein- 
Aether theory had been studied in [l[. Only the c\ term 
in plays an important role. To simplify the problem 
we choose C2 = C3 = 0. Since we are studying the weak 
gravity regime and the velocities of galaxies are much 
smaller than the speed of light, we assume that the time 
derivatives of the gravitational and vector fields are much 
smaller than their space derivatives. 

We expand the vector field and geometric metric 
around a fixed, Minkowski space background: 

9ap = VaP + h a p, 

A a = 5% +B a , 

where h a p and B a are taken to be small. We consider 
only terms that are linear in h and B, and choose Poisson 
gauge 



-2$ 



-2%*. 



(13) 



$ and ^ are scalar potentials. 

The first order term in equation (|12[) fixes 



B 


= 


as well as 




VtA* 
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S/tA* 


= a ( $ 


\J l A t 
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Vj* 


= d*B* 


K 


= M~ 2 



(14) 



We look for a solution such that diB^ = djB 1 , so that 
the first order contributions to the off-diagonal elements 
of the stress-energy tensor, T t i and Ty, disappear. The 
r/th component of the Einstein equations (i ^ j) then 
yields 



(15) 



This is consistent with General relativity. The 00th term 
of the Einstein equations is 

V 1 [(1 + y^'(K)) V* *] = 4*G N p. (16) 



This motivates us to identify 



M =l + |jr'(/C). 



(17) 



If diB' J = 0, then we recover the MOND form of grav- 
ity at small accelerations (i.e. at large distances). Here, 
we construct a solution that approaches the MOND be- 
havior at small acceleration even though diB^ ^ 0. We 
choose the functional form of 11 (really valid only for the 
acceleration due to a point mass) as in @ 



/' 



however, now 



((V$) 2 -(^') 2 ) 1 



1-p 



(18) 



(19) 



In principle, p is an arbitrary function of K, but we take 
it to be a constant. Note that 



(V$) 2 > {^f 



(20) 



so the presence of an inhomogeneous vector field estab- 
lishes a minimum acceleration. Moreover, the closer p is 
to 1, the greater the influence of the vector field. The 
usual MOND-type solution is recovered if 



($B J ') 2 =P(V3) S 



(21) 



We will take this condition as a starting point. Ulti- 
mately, one must demonstrate that the vector-field dy- 
namics leads to this relation, at least approximately. 
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METHOD 

Outside a single point mass the gravity in the deep 
MOND regime is 



V Gjvmao « 
V $ = r 



(22) 



The solution of the vector field for equation (f2"Tj) is then 



B = ± pGjfmao 



(23) 



Equation (|23)l is not well defined at r = 0; however, since 
we are not going to discuss the acceleration near the cen- 
ter, we ignore this divergence. 

Of more concern is that the vector field A a should 
approach (1,0,0,0), in other words B — ► 0, at spatial 
infinity. It seems that this requires that either the grav- 
ity deviates from the MOND limit (equation [22]) at very 
large r, or p(r) — ► at large r. At any rate, we shall 
assume that MOND applies on cluster scales, and that p 
remains constant on such scales, so that any deviations 
from these two approximations occurs on scales much 
larger than the size of a cluster. 

As has been shown in Q, perturbations of the vector 
field grow with the cosmological scale factor. Structure 
is formed as the vector field perturbations attract matter 
particles and thereby cause the gravitational potential 
to grow. The size and the mass of a galaxy or a clus- 
ter depends on the size and amplitude of the vector field 
perturbation. The fate of the seed vector perturbations, 
however, is unstudied. If these spatial terms of the vector 
field do not decay after the cluster structure is formed, 
then the space surrounding the cluster should still have 
a vector field that is different than the value that would 
be inferred by solving the vector field equation with the 
matter as a source. We assume that the vector field per- 
turbation extends over a region that is larger than the 
size of our simulation, about 8Mpc. 

For two point particles with mass mi and m.2, the grav- 
itational acceleration in the MOND regime superposes 
as: 



(v*) 2 ~ (V*i) 2 + (v$ 2 ) 2 = 



(24) 

Note that this is very different than the usual Newtonian 
superposition of V*£- To satisfy this relation we assume 
that the vector fields adding according to: 

(d t BJ) 2 = (diB{) 2 + {diB{f, 

where B\ and Bi are the vector fields associated with 
(but not necessarily sourced by) m\ and m; . 

There are several numerical methods 17, li|, 2(| that 
can be used to calculate the gravitational acceleration 
in the environments of a mass distribution according to 



MOND. To simplify this problem, we choose an axially 
symmetric mass distribution. The method in 17[ is more 
suitable and easier for this case. The code solves the 
following equations: 



Here 



y -U = AttGnP 
V x vU = 0. 



U = v * 



(25) 
(26) 



(27) 
(28) 



serve as definitions of U and v respecitively. Here, /i is 
defined by equations (fl"8|) and (Tj"9|) . On a cylindrical 



grid, defined by {(rj, Zj)} (where i indexes the radial co- 
ordinate, r, and j indexes the axial coordinate z, and 
the symmetry allows us to ignore the azimuthal coordi- 
nate), these differential equations can be replaced by the 
following difference equations: 



(U, 



-i u i+i,j - r i U i,j) =- 4 ' Y '" 



-2it5zj(ri. 



(29) 



= 0. 



(30) 



Here U[j is the r component of U at 7"j, Zj\ similarly £7f 



is the z component. The r)k are discrete representations 
of fxdl, eg. rji = (i/j_i5rj_i + Vij8ri)j2. For details of 
the calculational method, please refer to [l7T |. 

Since our goal is to investigate whether the weak lens- 
ing maps of the bullet cluster can be reproduced in the 
MOND-like limit of GEA, we must calculate those maps. 
The lensing convergence is 



DdsDd 
c 2 D s 



V 2 <S>dl 



(31) 



Here, Dd, Dd s and D s are the angular distances between 
the observer and the lens, the lens and the source, and the 
observer and the source, respectively, i is the distance 
along the light ray from the observer. Once we have 
solved equations (|29p and (130)) , we can straightforwardly 
calculate k. 



MODEL AND RESULT 

Our code is able to simulate only axially symmetric 
mass distributions. The mass distribution of the bullet 
cluster, fortunately, is very close to axially symmetric. 
The angle between the direction of the motion of the 
bullet cluster and the the plane of the sky is about 8°. 
We, therefore, assume that the bullet cluster is axially 
symmetric about an axis that is parallel to the plane of 
the sky. 
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FIG. 2: The covergence map changing with different p: — 50kpc for the vector fields, p = 0.00, p — 0.50, p = 0.9, and 
p = 0.99 are corresponding to left up, right up, left bottom, and right bottom. The interval between each contour is 0.15. 
Left up figure shows the result of MOND. It does not show a gravitational lensing center at — 350kpc and 350fcpc. As the p 
increases, concentracted contours appears. 



For the baryonic mass distribution of the cluster, we 
take a King profile [28[ : 



p(r) 



Po 



i + (£) 5 



(32) 



We choose a core radius, r c = 150kpc. The central den- 
sity, po, and the positions of mass centers are chosen 
acording to table Q] as extracted from 21] . Figure Q] We 
choose l.hMpc to be the cutoff of the mass distribution. 
Dd £ ) Dd = 540-Mpc, which is appropriate for the lensed 
galaxies in the background of the bullet cluster, where the 
source of the gravitational lensing is at Zi — 0.296, and 
the a mean redshift of the lensed galaxies is Zb g = 0.85 

In addition to the baryonic clusters, we add vector field 
configurations at position A and D. Since the only cou- 
pling constant can appear in our model is Gn and 
ao both of them is very small. Therefore, we as- 
sume the interaction between vector field is small. 



It can pass through each other unimpeded. The 

mass term in equation (f2"2"|) for the vector field at A is cho- 
sen to be the mass of component A plus that of compo- 
nent B. Similarly, the mass term in equation (|22[) for the 
vector field at D is chosen to be the mass of component 
C plus that of component D. We assume the mass also 
follows equation (|32|) . But the r^s of the vector fields are 
different from those of the baryonic mass distributions. 

Figure [5] shows the convergence map changing as p 
is varied. Here, we choose = bOkpc for the vector 
fields. For p = 0, the solution is the same as MOND 
There are two lensing centers that appear in the map, 
but they are located approximately halfway between A 
and B, and halfway between C and D. Thus, to repro- 
duce the lensing map of the bullet cluster would require 
vector-concentrations that are very far outside the bary- 
onic mass distributions. 

However, as p is increased, the centers of the conver- 
gence contours move away from the baryonic centers, and 
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FIG. 3: The covergence map changing with different for the vector fields: p = 0.99 for the model, — 200kpc, — lOOfcpc, 
th = 25kpc, and ru = Wkpc are corresponding to left up, right up, left bottom, and right bottom. The interval between each 
contour is 0.15. As the decreases, concentracted contours appears near — 350fcpc and 350kpc. 



TABLE I: The mass and position of baryon matter 



component 


type 


position 
(kpc) 


Plasma mass 
within 100/spc 


A 


lensing 
center 


-350 


6.0 xlO 12 M 


B 


X-ray 
center 


-150 


6.8 xl0 12 M© 


C 


X-ray 
center 


200 


3.3 xl0 12 M Q 


D 


lensing 
center 


350 


5.9 xlO 12 M 



k is much larger than in MOND. This means that for p 
near 1, the vector fields play a role very much like dark 
matter, and are likely to augment the lensing in very 
much the same way as do dark matter halos. 

Figure [3] shows the convergence map changing as 



for the vector fields is varied. To emphasize the effect of 
the vector field, we choose p = 0.99. Moving from panel 
a (top left) through panel d (bottom right), changes 
from 200fcpc to lOkpc. For large r^, k is less than in 
MOND. This is because (from equation [T9| the effective 
ao is 



a e ff = a \J\ - p. (33) 

The MOND acceleration regime applies only when a < 
a e ff. If the scale height of the vector field (ru) is very 
large, then a < a e ff only at very large distances; at 
smaller distances, the gravity becomes Newtonian, k is 
then much less than it would have been under MOND. As 
r/j is reduced, a gravitational center appears and moves 
to the centers of the vector fields and the n also becomes 
bigger. Clearly, Figure [5] and [3] suggest that the space 
components of vector fields can replace dark matter. 
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CONCLUSION 

In this work, we analyzed the effects of gravitational 
lcnsing caused by the vector field in the Generalized 
Einstein- Aether theory. The results show that the vec- 
tor field could cause the dissociation of mass center and 
gravity center observed in the bullet cluster. This would 
not be the vector field sourced by the matter, but rather 
the vector field that seeded the growth of the matter 
perturbation. We have not followed the evolution of this 
primordial seed, but rather established the approximate 
properties that it would have to explain the observations. 

The problem with the vector field is that a particular 
configuration of the vector field distribution is needed. 
This is the same problem as the need for a particular mass 
distribution of the dark matter. Future work will address 
comparing the full weak lensing data instead of the con- 
vergence map only. We note that, since TeVeS also in- 
cludes the same type of vector field, it is reasonable to 
believe that the gravitational acceleration of TeVeS could 
be dissociated from mass centers by the vector fields. 

The convergence map alone is not enough to distin- 
guish this type of vector field from dark matter. We 
need a more detailed analysis to distinguish the two mod- 
els, eg. by looking at the lensing shear we may yet find 
that modified gravity is distinguishable from dark mat- 
ter. Finally we note that this independent concentration 
of vector is very much like a dark matter halo - but a 
halo of vector field rather than of particles. Thus MOND 
started as an effort to build a theory in which what you 
see (baryons) is what you get (i.e. sources all gravity). 
One way or another, even in modified gravity theories, 
that seems not to be true. 

It is a pleasure to thank Constantinos Skordis, Irit 
Maor, Pedro Ferreira and Thomas Zlosnik for very useful 
discussions. We thank the authors of 2l| for permission 
to reproduce figure [TJ This work has been supported by 
grants from the US DOE, NASA and HEPCOS group at 
SUNY. 
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We study the lensing properties of an asymmetric mass distribution and vector field in General- 
ized Einstein-Aether (GEA) theory. As vector field fluctuations are responsible in GEA for seeding 
baryonic structure formation, vector field concentrations can exist independently of baryonic mat- 
ter. Such concentrations would not be expected to be tied to baryonic matter except gravitationally, 
and so, like dark matter halos, would become separated from baryonic matter in interacting systems 
such as the Bullet Cluster. These vector field concentrations cause metric deviations that affect 
weak lensing. Therefore, the distribution of weak lensing deviates from that which would be in- 
ferred from the luminous mass distribution, in a way that numerical calculations demonstrate can 
be consistent with observations. This suggests that MOND-inspired theories can reproduce weak 
lensing observations, but makes clear the price: the existence of a coherent large-scale fluctuation 
of a field(s) weakly tied to the baryonic matter, not completely dissimilar to a dark matter halo. 
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INTRODUCTION 

It is well known that General Relativity can not explain 
the full dynamics of our universe and the structures in it 
sourced solely by visible mass. Two kinds of exotic energy 
density, dark matter and dark energy, must be introduced 
to explain respectively the dynamics of structures and the 
accelerating expansion of the universe. 

An alternative to adding new sources of gravity is to 
alter the response of gravity to the known matter sources. 
Modified theories of gravity replacing dark matter as an 
explanation of flat galaxy rotation curves were first pro- 
posed by Milgrom [fl, 0] and by Milgrom and Bekenstein 

@] , and came to be known as Modified Newtonian Dy- 
namics a.k.a. MOND (for a review of MOND see 0). In 
the contemporary incarnation of this phenomenological 
theory, Poisson's equation for the Newtonian potential 
(j), is altered by the inclusion of a scalar function fi: 



V-(M(rW>|)V£) =4*G N p. 



(1) 



The scalar function is chosen to enhance gravitational 
acceleration at large distance, and thus replace galactic 
dark matter. It turns out that galactic dynamics be ex- 
plained using very simple choices of /x, such as 



m(IW|) 



|V<^|/a 



(2) 



Introducing only one new parameter, a universal acceler- 
ation, do, one can fit almost all galactic rotation curves 
assuming very reasonable mass-to-light ratios consistent 
with prior expectations. 

Despite its successes, most scientists did not regard 
MOND as a viable theory. In part this was due to 
difficulties explaining clusters and possibly some small 



scale structures, but more particularly it was due to 
the absence of a consistent covariant action formula- 
tion of the theory. More recently, these theoretical hur- 
dles were overcome by Bekenstein who proposed 
Tensor- Vector-Scalar gravity (TeVeS). Zlosnik et al.^ 
then showed that TeVeS is part of a larger class that 
they call Generalized Einstein- Aether (GEA) theories. In 
these theories, a new field or fields modify the response 
of the metric to the presence of matter. Most generic is 
the existence of a vector field with a non-zero time-like 
background value. Hence the name Einstein-Aether. 

It has been claimed that the discovery of the merging 
galaxy cluster 1E0657-56 the bullet cluster, rules 
out these modified gravity theories and proves the exis- 
tence of dark matter [2]]]. X-ray images of this cluster 
(left hand side of figure [T]) reveal a bullet-like subcluster 
just exiting the core of the main cluster. Weak lensing 
mass contours (right hand side of Figure [T]) show that 
the lensing centers are not at the luminous centers of 
these bright clusters, rather they are (approximately) on 
the line through these X-ray centers but outside them 
(cf. Figure [1]). The dark matter interpretation is that 
the lensing centers (which are less luminous than the 
X-ray centers) are caused by the presence of collision- 
less cold dark matter (CDM). In the aftermath of the 
cluster-cluster collision, the CDM halos have presumably 
continued along un-impeded, while the hot baryonic gas 
has been slowed. 

It is true that MOND is unlikely to explain the sepa- 
ration of the gravity centers from the mass centers. But 
MOND was never to be regarded as a complete the- 
ory, and certainly not as a theory of time-dependent 
mass distributions. Previous authors 2H, 2fi| have taken 
the MOND-derived Newtonian potential, inserted it into 
TeVeS to show that TeVeS too does not explain the ob- 
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FIG. 1: The bullet cluster [2l|: Shown above in the left panel is a color image from the Magellan images of the merging cluster 
1E0657— 558, with the white bar indicating 200 kpc at the distance of the cluster. In the right panel is a 500 ks Chandra image 
of the cluster. Shown in green contours in both panels are the weak lensing k reconstruction with the outer contour level at 
ft = 0.16 and increasing in steps of 0.07. The white contours show the errors on the positions of the k peaks and correspond 
to 68.3%, 95.5%, and 99.7% confidence levels. 



served weak lensing map without an additional dark mass 
component. This has led the Bullet cluster to be widely 
regarded as direct evidence for dark matter and against 
modified gravity. 

But are relativistic versions of MOND truly ruled out 
by the bullet cluster? The new dynamical degrees of 
freedom of these theories are known to solve other clear 
problems of vanilla-MOND. For example, it was shown 
[lfj that perturbations glows too slowly in a MOND 
universe for the observed large scale structures to have 
emerged from the fluctuations present at the time of re- 
combination. However, it has been shown in both TeVeS 



111 [12| and GEA growing modes of the new dynamical 
fields introduced in these modified gravity theories suc- 
cessfully seed the growth of baryonic structures after the 
end of silk damping. This implies that observed bary- 
onic structures - galaxies, clusters, etc. - evolved from 
non-baryonic seeds which were not sourced by baryons. 
If these seeds persist, then they might be expected to 
behave in very much the same way as a dark matter 
halo when the baryonic structures with which they are 
associated collide - travelling on, in a straight line unim- 
peded. These non-baryonic structures may then be the 
non-luminous weak lensing centers observed in the bullet 
cluster. Of course, this would mean that even in modi- 
fied gravity theories baryonic structures are surrounded 
by non-baryonic halos - but halos of weakly interacting 
classical fields rather than halos of weakly interacting 
non-relativistic particles. 

For our calculations, we will focus on the Generalized 
Einstein- Aether theory, a generalized form of " Einstern- 
Aether" theory [3) El • I n this theory, general relativity 
is modified only by the addition of a dynamical vector 
field A a , with a time- like classical value A a = (1, 0, 0, 0). 
We modify the conventional ansatz for A a in the presence 
of a matter source to allow for first-order perturbations in 



the spacelike, rather than just the time-like components. 
The spacelike component is not sourced by the matter, 
but could emerge from the intrinsic vector growing mode. 
The vector field enhances weak lensing by different fac- 
tors at different locations, with a resulting possible disso- 
ciation of the luminous and gravitational-lensing centers. 
MOG is another possible modified gravity theory htat 
mijg ht accomodate the observations of the bullet cluster 

We organize the paper as follows. We first review 
briefly the Generalized Einstein- Aether theory and cal- 
culate the solution for the vector field to first order. We 
then discuss appropriate boundary conditions for the bul- 
let cluster, and perform suitable numerical calculations 
of the resulting metric and convergence. 



THEORY 

The action of the Generalized Einstein- Aether theory 
can be written in the form [l[ 



S = 



R 



16ttG 



A' 



£(9, A) 



S 



M ■ 



(3) 



Here, g is the usual metric (with signature (-,+,+,+)), R 
the Ricci scalar of that metric, and Sm the matter action. 
A is a new dynamical vector field. As shown in [H, @], 
we can obtain a MOND-type limit in the appropriate 
low-acceleration regime by choosing 



where 

K, 



M 2 



16ttG 



1 



JV 



16ttG 



N 



-X(A a A a + 1), (4) 



K 



M- 2 K. al3 V a A~<V B A a and 



a/3 



75 



(5) 
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The Ci are dimensionless constants, whereas M has the 
dimension of mass. A is a non-dynamical Lagrange- 
multiplier field with dimensions of mass-squared, that 
enforces that A a is unit-timelike. The gravitational field 
equations obtained by varying g a @ @, [jj) are 



G a p — T a p + 87rGArT a/3 



(6) 



where the stress-energy tensor for the vector field is given 
by 



1 



-T'Y< 



(a/3) 



-g a pM 2 T + XA a A 



(7) 



where 



T = die and 

(/C Q/3 ^ 



(8) 



Brackets around indices denote symmetrization. Y a p is 
the functional derivative 



Y 



a/3 



= V.A^A^^ 



7 ^ 6g a P 
= -ciKVMiVA^-iVMiVpA")} 

The equations obtained by varying A@ are 

V a {T l J a l3 ) + T' V p = 2XAp, 

with 



(9) 



(10) 



SAP 



For the choice of K. given by equation [H yp = 0. 
Variations of A fix 



A^A, 



-1. 



(11) 



(12) 



The MOND type solutions of the Generalized Einstein- 
Aether theory had been studied in [l[. Only the c\ term 
in plays an important role. To simplify the problem 
we choose C2 = C3 = 0. Since we are studying the weak 
gravity regime and the velocities of galaxies are much 
smaller than the speed of light, we assume that the time 
derivatives of the gravitational and vector fields are much 
smaller than their space derivatives. 

We expand the vector field and geometric metric 
around a fixed, Minkowski space background: 

9ap = VaP + h a p, 

A a = 5% +B a , 

where h a p and B a are taken to be small. We consider 
only terms that are linear in h and B, and choose Poisson 
gauge 



-2$ 



-2%*. 



(13) 



$ and ^ are scalar potentials. 

The first order term in equation (|12[) fixes 



B 


= 


as well as 




VtA* 


= 


S/tA* 


= a ( $ 


\J l A t 


= 


Vj* 


= d*B* 


K 


= M~ 2 



(14) 



We look for a solution such that diB^ = djB 1 , so that 
the first order contributions to the off-diagonal elements 
of the stress-energy tensor, T t i and Ty, disappear. The 
r/th component of the Einstein equations (i ^ j) then 
yields 



(15) 



This is consistent with General relativity. The 00th term 
of the Einstein equations is 

V 1 [(1 + y^'(K)) V* *] = 4*G N p. (16) 



This motivates us to identify 



M =l + |jr'(/C). 



(17) 



If diB' J = 0, then we recover the MOND form of grav- 
ity at small accelerations (i.e. at large distances). Here, 
we construct a solution that approaches the MOND be- 
havior at small acceleration even though diB^ ^ 0. We 
choose the functional form of 11 (really valid only for the 
acceleration due to a point mass) as in @ 



/' 



however, now 



((V$) 2 -(^') 2 ) 1 



1-p 



(18) 



(19) 



In principle, p is an arbitrary function of K, but we take 
it to be a constant. Note that 



(V$) 2 > {^f 



(20) 



so the presence of an inhomogeneous vector field estab- 
lishes a minimum acceleration. Moreover, the closer p is 
to 1, the greater the influence of the vector field. The 
usual MOND-type solution is recovered if 



($B J ') 2 =P(V3) S 



(21) 



We will take this condition as a starting point. Ulti- 
mately, one must demonstrate that the vector-field dy- 
namics leads to this relation, at least approximately. 
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METHOD 

Outside a single point mass the gravity in the deep 
MOND regime is 



V Gjvmao « 
V $ = r 



(22) 



The solution of the vector field for equation (f2"Tj) is then 



B = ± pGjfmao 



(23) 



Equation (|23)l is not well defined at r = 0; however, since 
we are not going to discuss the acceleration near the cen- 
ter, we ignore this divergence. 

Of more concern is that the vector field A a should 
approach (1,0,0,0), in other words B — ► 0, at spatial 
infinity. It seems that this requires that either the grav- 
ity deviates from the MOND limit (equation [22]) at very 
large r, or p(r) — ► at large r. At any rate, we shall 
assume that MOND applies on cluster scales, and that p 
remains constant on such scales, so that any deviations 
from these two approximations occurs on scales much 
larger than the size of a cluster. 

As has been shown in Q, perturbations of the vector 
field grow with the cosmological scale factor. Structure 
is formed as the vector field perturbations attract matter 
particles and thereby cause the gravitational potential 
to grow. The size and the mass of a galaxy or a clus- 
ter depends on the size and amplitude of the vector field 
perturbation. The fate of the seed vector perturbations, 
however, is unstudied. If these spatial terms of the vector 
field do not decay after the cluster structure is formed, 
then the space surrounding the cluster should still have 
a vector field that is different than the value that would 
be inferred by solving the vector field equation with the 
matter as a source. We assume that the vector field per- 
turbation extends over a region that is larger than the 
size of our simulation, about 8Mpc. 

For two point particles with mass mi and m.2, the grav- 
itational acceleration in the MOND regime superposes 
as: 



(v*) 2 ~ (V*i) 2 + (v$ 2 ) 2 = 



(24) 

Note that this is very different than the usual Newtonian 
superposition of V*£- To satisfy this relation we assume 
that the vector fields adding according to: 

(d t BJ) 2 = (diB{) 2 + {diB{f, 

where B\ and Bi are the vector fields associated with 
(but not necessarily sourced by) m\ and m; . 

There are several numerical methods 17, li|, 2(| that 
can be used to calculate the gravitational acceleration 
in the environments of a mass distribution according to 



MOND. To simplify this problem, we choose an axially 
symmetric mass distribution. The method in 17[ is more 
suitable and easier for this case. The code solves the 
following equations: 



Here 



y -U = AttGnP 
V x vU = 0. 



U = v * 



(25) 
(26) 



(27) 
(28) 



serve as definitions of U and v respecitively. Here, /i is 
defined by equations (fl"8|) and (Tj"9|) . On a cylindrical 



grid, defined by {(rj, Zj)} (where i indexes the radial co- 
ordinate, r, and j indexes the axial coordinate z, and 
the symmetry allows us to ignore the azimuthal coordi- 
nate), these differential equations can be replaced by the 
following difference equations: 



(U, 



-i u i+i,j - r i U i,j) =- 4 ' Y '" 



-2it5zj(ri. 



(29) 



= 0. 



(30) 



Here U[j is the r component of U at 7"j, Zj\ similarly £7f 



is the z component. The r)k are discrete representations 
of fxdl, eg. rji = (i/j_i5rj_i + Vij8ri)j2. For details of 
the calculational method, please refer to [l7T |. 

Since our goal is to investigate whether the weak lens- 
ing maps of the bullet cluster can be reproduced in the 
MOND-like limit of GEA, we must calculate those maps. 
The lensing convergence is 



DdsDd 
c 2 D s 



V 2 <S>dl 



(31) 



Here, Dd, Dd s and D s are the angular distances between 
the observer and the lens, the lens and the source, and the 
observer and the source, respectively, i is the distance 
along the light ray from the observer. Once we have 
solved equations (|29p and (130)) , we can straightforwardly 
calculate k. 



MODEL AND RESULT 

Our code is able to simulate only axially symmetric 
mass distributions. The mass distribution of the bullet 
cluster, fortunately, is very close to axially symmetric. 
The angle between the direction of the motion of the 
bullet cluster and the the plane of the sky is about 8°. 
We, therefore, assume that the bullet cluster is axially 
symmetric about an axis that is parallel to the plane of 
the sky. 
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FIG. 2: The covergence map changing with different p: — 50kpc for the vector fields, p = 0.00, p — 0.50, p = 0.9, and 
p = 0.99 are corresponding to left up, right up, left bottom, and right bottom. The interval between each contour is 0.15. 
Left up figure shows the result of MOND. It does not show a gravitational lensing center at — 350kpc and 350fcpc. As the p 
increases, concentracted contours appears. 



For the baryonic mass distribution of the cluster, we 
take a King profile [28[ : 



p(r) 



Po 



i + (£) 5 



(32) 



We choose a core radius, r c = 150kpc. The central den- 
sity, po, and the positions of mass centers are chosen 
acording to table Q] as extracted from 21] . Figure Q] We 
choose l.hMpc to be the cutoff of the mass distribution. 
Dd £ ) Dd = 540-Mpc, which is appropriate for the lensed 
galaxies in the background of the bullet cluster, where the 
source of the gravitational lensing is at Zi — 0.296, and 
the a mean redshift of the lensed galaxies is Zb g = 0.85 

In addition to the baryonic clusters, we add vector field 
configurations at position A and D. The only coupling 
constants in the model are the usual gravitational cou- 
pling Gn, and the very small acceleration aQ. We there- 
fore assume that the interaction between the vector field 



"halos" associated with each cluster are sufficiently small 
that the two vector-halos can pass through each other 
unimpeded, much as dark matter halos are assumed to. 
The mass term in equation (|22p for the vector field at 
A is chosen to be the mass of component A plus that 
of component B. Similarly, the mass term in equation 
(f2"2"]l for the vector field at D is chosen to be the mass 
of component C plus that of component D. We assume 
the mass also follows equation ([3"2"]l . But the r^s of the 
vector fields are different from those of the baryonic mass 
distributions. 

Figure [5] shows the convergence map changing as p 
is varied. Here, we choose = bOkpc for the vector 
fields. For p = 0, the solution is the same as MOND 
There are two lensing centers that appear in the map, 
but they are located approximately halfway between A 
and B, and halfway between C and D. Thus, to repro- 
duce the lensing map of the bullet cluster would require 
vector-concentrations that are very far outside the bary- 
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FIG. 3: The convergence map changing with different for the vector fields: p = 0.99 for the model, — 200kpc, — lOOfcpc, 
th = 25kpc, and ru = Wkpc are corresponding to left up, right up, left bottom, and right bottom. The interval between each 
contour is 0.15. As the decreases, concentracted contours appears near — 350fcpc and 350£;pc. 



TABLE I: The mass and position of baryon matter 



component 


type 


position 
(kpc) 


Plasma mass 
within 100/spc 


A 


lensing 
center 


-350 


6.0 xlO 12 M 


B 


X-ray 
center 


-150 


6.8 xl0 12 M© 


C 


X-ray 
center 


200 


3.3 xl0 12 M Q 


D 


lensing 
center 


350 


5.9 xlO 12 M 



onic mass distributions. 

However, as p is increased, the centers of the conver- 
gence contours move away from the baryonic centers, and 
k is much larger than in MOND. This means that for p 
near 1, the vector fields play a role very much like dark 



matter, and are likely to augment the lensing in very 
much the same way as do dark matter halos. 

Figure [3] shows the convergence map changing as Th 
for the vector fields is varied. To emphasize the effect of 
the vector field, we choose p = 0.99. Moving from panel 
a (top left) through panel d (bottom right), changes 
from 200fcpc to lOkpc. For large r^, k is less than in 
MOND. This is because (from equation [T9f the effective 
ao is 

a e ff = a Q ^l -p- (33) 

The MOND acceleration regime applies only when a < 
a e ff. If the scale height of the vector field (ru) is very 
large, then a < a e ff only at very large distances; at 
smaller distances, the gravity becomes Newtonian, k is 
then much less than it would have been under MOND. As 
r/j is reduced, a gravitational center appears and moves 
to the centers of the vector fields and the n also becomes 
bigger. Clearly, Figure [5] and [3] suggest that the space 
components of vector fields can replace dark matter. 
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CONCLUSION 

In this work, we analyzed the effects of gravitational 
lcnsing caused by the vector field in the Generalized 
Einstein- Aether theory. The results show that the vec- 
tor field could cause the dissociation of mass center and 
gravity center observed in the bullet cluster. This would 
not be the vector field sourced by the matter, but rather 
the vector field that seeded the growth of the matter 
perturbation. We have not followed the evolution of this 
primordial seed, but rather established the approximate 
properties that it would have to explain the observations. 

The problem with the vector field is that a particular 
configuration of the vector field distribution is needed. 
This is the same problem as the need for a particular mass 
distribution of the dark matter. Future work will address 
comparing the full weak lensing data instead of the con- 
vergence map only. We note that, since TeVeS also in- 
cludes the same type of vector field, it is reasonable to 
believe that the gravitational acceleration of TeVeS could 
be dissociated from mass centers by the vector fields. 

The convergence map alone is not enough to distin- 
guish this type of vector field from dark matter. We 
need a more detailed analysis to distinguish the two mod- 
els, eg. by looking at the lensing shear we may yet find 
that modified gravity is distinguishable from dark mat- 
ter. Finally we note that this independent concentration 
of vector is very much like a dark matter halo - but a 
halo of vector field rather than of particles. Thus MOND 
started as an effort to build a theory in which what you 
see (baryons) is what you get (i.e. sources all gravity). 
One way or another, even in modified gravity theories, 
that seems not to be true. 

It is a pleasure to thank Constantinos Skordis, Irit 
Maor, Pedro Ferreira and Thomas Zlosnik for very useful 
discussions. We thank the authors of 2l| for permission 
to reproduce figure [TJ This work has been supported by 
grants from the US DOE, NASA and HEPCOS group at 
SUNY. 
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We study the lensing properties of an asymmetric mass distribution and vector field in General- 
ized Einstein-Aether (GEA) theory. As vector field fluctuations are responsible in GEA for seeding 
baryonic structure formation, vector field concentrations can exist independently of baryonic mat- 
ter. Such concentrations would not be expected to be tied to baryonic matter except gravitationally, 
and so, like dark matter halos, would become separated from baryonic matter in interacting systems 
such as the Bullet Cluster. These vector field concentrations cause metric deviations that affect 
weak lensing. Therefore, the distribution of weak lensing deviates from that which would be in- 
ferred from the luminous mass distribution, in a way that numerical calculations demonstrate can 
be consistent with observations. This suggests that MOND-inspired theories can reproduce weak 
lensing observations, but makes clear the price: the existence of a coherent large-scale fluctuation 
of a field(s) weakly tied to the baryonic matter, not completely dissimilar to a dark matter halo. 
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INTRODUCTION 

It is well known that General Relativity can not explain 
the full dynamics of our universe and the structures in it 
sourced solely by visible mass. Two kinds of exotic energy 
density, dark matter and dark energy, must be introduced 
to explain respectively the dynamics of structures and the 
accelerating expansion of the universe. 

An alternative to adding new sources of gravity is to 
alter the response of gravity to the known matter sources. 
Modified theories of gravity replacing dark matter as an 
explanation of flat galaxy rotation curves were first pro- 
posed by Milgrom [fl, 0] and by Milgrom and Bekenstein 

@] , and came to be known as Modified Newtonian Dy- 
namics a.k.a. MOND (for a review of MOND see 0). In 
the contemporary incarnation of this phenomenological 
theory, Poisson's equation for the Newtonian potential 
(j), is altered by the inclusion of a scalar function fi: 



V-(M(rW>|)V£) =4*G N p. 



(1) 



The scalar function is chosen to enhance gravitational 
acceleration at large distance, and thus replace galactic 
dark matter. It turns out that galactic dynamics be ex- 
plained using very simple choices of /x, such as 



m(IW|) 



|V<^|/a 



(2) 



Introducing only one new parameter, a universal acceler- 
ation, do, one can fit almost all galactic rotation curves 
assuming very reasonable mass-to-light ratios consistent 
with prior expectations. 

Despite its successes, most scientists did not regard 
MOND as a viable theory. In part this was due to 
difficulties explaining clusters and possibly some small 



scale structures, but more particularly it was due to 
the absence of a consistent covariant action formula- 
tion of the theory. More recently, these theoretical hur- 
dles were overcome by Bekenstein who proposed 
Tensor- Vector-Scalar gravity (TeVeS). Zlosnik et al.^ 
then showed that TeVeS is part of a larger class that 
they call Generalized Einstein- Aether (GEA) theories. In 
these theories, a new field or fields modify the response 
of the metric to the presence of matter. Most generic is 
the existence of a vector field with a non-zero time-like 
background value. Hence the name Einstein-Aether. 

It has been claimed that the discovery of the merging 
galaxy cluster 1E0657-56 the bullet cluster, rules 
out these modified gravity theories and proves the exis- 
tence of dark matter [2]]]. X-ray images of this cluster 
(left hand side of figure [T]) reveal a bullet-like subcluster 
just exiting the core of the main cluster. Weak lensing 
mass contours (right hand side of Figure [T]) show that 
the lensing centers are not at the luminous centers of 
these bright clusters, rather they are (approximately) on 
the line through these X-ray centers but outside them 
(cf. Figure [1]). The dark matter interpretation is that 
the lensing centers (which are less luminous than the 
X-ray centers) are caused by the presence of collision- 
less cold dark matter (CDM). In the aftermath of the 
cluster-cluster collision, the CDM halos have presumably 
continued along un-impeded, while the hot baryonic gas 
has been slowed. 

It is true that MOND is unlikely to explain the sepa- 
ration of the gravity centers from the mass centers. But 
MOND was never to be regarded as a complete the- 
ory, and certainly not as a theory of time-dependent 
mass distributions. Previous authors 2H, 2fi| have taken 
the MOND-derived Newtonian potential, inserted it into 
TeVeS to show that TeVeS too does not explain the ob- 
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FIG. 1: The bullet cluster [2l|: Shown above in the left panel is a color image from the Magellan images of the merging cluster 
1E0657— 558, with the white bar indicating 200 kpc at the distance of the cluster. In the right panel is a 500 ks Chandra image 
of the cluster. Shown in green contours in both panels are the weak lensing k reconstruction with the outer contour level at 
ft = 0.16 and increasing in steps of 0.07. The white contours show the errors on the positions of the k peaks and correspond 
to 68.3%, 95.5%, and 99.7% confidence levels. 



served weak lensing map without an additional dark mass 
component. This has led the Bullet cluster to be widely 
regarded as direct evidence for dark matter and against 
modified gravity. 

But are relativistic versions of MOND truly ruled out 
by the bullet cluster? The new dynamical degrees of 
freedom of these theories are known to solve other clear 
problems of vanilla-MOND. For example, it was shown 
[lfj that perturbations glows too slowly in a MOND 
universe for the observed large scale structures to have 
emerged from the fluctuations present at the time of re- 
combination. However, it has been shown in both TeVeS 



111 [12| and GEA growing modes of the new dynamical 
fields introduced in these modified gravity theories suc- 
cessfully seed the growth of baryonic structures after the 
end of silk damping. This implies that observed bary- 
onic structures - galaxies, clusters, etc. - evolved from 
non-baryonic seeds which were not sourced by baryons. 
If these seeds persist, then they might be expected to 
behave in very much the same way as a dark matter 
halo when the baryonic structures with which they are 
associated collide - travelling on, in a straight line unim- 
peded. These non-baryonic structures may then be the 
non-luminous weak lensing centers observed in the bullet 
cluster. Of course, this would mean that even in modi- 
fied gravity theories baryonic structures are surrounded 
by non-baryonic halos - but halos of weakly interacting 
classical fields rather than halos of weakly interacting 
non-relativistic particles. 

For our calculations, we will focus on the Generalized 
Einstein- Aether theory, a generalized form of " Einstern- 
Aether" theory [3) El • I n this theory, general relativity 
is modified only by the addition of a dynamical vector 
field A a , with a time- like classical value A a = (1, 0, 0, 0). 
We modify the conventional ansatz for A a in the presence 
of a matter source to allow for first-order perturbations in 



the spacelike, rather than just the time-like components. 
The spacelike component is not sourced by the matter, 
but could emerge from the intrinsic vector growing mode. 
The vector field enhances weak lensing by different fac- 
tors at different locations, with a resulting possible disso- 
ciation of the luminous and gravitational-lensing centers. 
MOG is another possible modified gravity theory htat 
mijg ht accomodate the observations of the bullet cluster 

We organize the paper as follows. We first review 
briefly the Generalized Einstein- Aether theory and cal- 
culate the solution for the vector field to first order. We 
then discuss appropriate boundary conditions for the bul- 
let cluster, and perform suitable numerical calculations 
of the resulting metric and convergence. 



THEORY 

The action of the Generalized Einstein- Aether theory 
can be written in the form [l[ 



S = 



R 



16ttG 



A' 



£(9, A) 



S 



M ■ 



(3) 



Here, g is the usual metric (with signature (-,+,+,+)), R 
the Ricci scalar of that metric, and Sm the matter action. 
A is a new dynamical vector field. As shown in [H, @], 
we can obtain a MOND-type limit in the appropriate 
low-acceleration regime by choosing 



where 

K, 



M 2 



16ttG 



1 



JV 



16ttG 



N 



-X(A a A a + 1), (4) 



K 



M- 2 K. al3 V a A~<V B A a and 



a/3 



75 



(5) 
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The Ci are dimensionless constants, whereas M has the 
dimension of mass. A is a non-dynamical Lagrange- 
multiplier field with dimensions of mass-squared, that 
enforces that A a is unit-timelike. The gravitational field 
equations obtained by varying g a @ @, [jj) are 



G a p — T a p + 87rGArT a/3 



(6) 



where the stress-energy tensor for the vector field is given 
by 



1 



-T'Y< 



(a/3) 



-g a pM 2 T + XA a A 



(7) 



where 



T = die and 

(/C Q/3 ^ 



(8) 



Brackets around indices denote symmetrization. Y a p is 
the functional derivative 



Y 



a/3 



= V.A^A^^ 



7 ^ 6g a P 
= -ciKVMiVA^-iVMiVpA")} 

The equations obtained by varying A@ are 

V a {T l J a l3 ) + T' V p = 2XAp, 

with 



(9) 



(10) 



SAP 



For the choice of K. given by equation [H yp = 0. 
Variations of A fix 



A^A, 



-1. 



(11) 



(12) 



The MOND type solutions of the Generalized Einstein- 
Aether theory had been studied in [l[. Only the c\ term 
in plays an important role. To simplify the problem 
we choose C2 = C3 = 0. Since we are studying the weak 
gravity regime and the velocities of galaxies are much 
smaller than the speed of light, we assume that the time 
derivatives of the gravitational and vector fields are much 
smaller than their space derivatives. 

We expand the vector field and geometric metric 
around a fixed, Minkowski space background: 

9ap = VaP + h a p, 

A a = 5% +B a , 

where h a p and B a are taken to be small. We consider 
only terms that are linear in h and B, and choose Poisson 
gauge 



-2$ 



-2%*. 



(13) 



$ and ^ are scalar potentials. 

The first order term in equation (|12[) fixes 



B 


= 


as well as 




VtA* 


= 


S/tA* 


= a ( $ 


\J l A t 


= 


Vj* 


= d*B* 


K 


= M~ 2 



(14) 



We look for a solution such that diB^ = djB 1 , so that 
the first order contributions to the off-diagonal elements 
of the stress-energy tensor, T t i and Ty, disappear. The 
r/th component of the Einstein equations (i ^ j) then 
yields 



(15) 



This is consistent with General relativity. The 00th term 
of the Einstein equations is 

V 1 [(1 + y^'(K)) V* *] = 4*G N p. (16) 



This motivates us to identify 



M =l + |jr'(/C). 



(17) 



If diB' J = 0, then we recover the MOND form of grav- 
ity at small accelerations (i.e. at large distances). Here, 
we construct a solution that approaches the MOND be- 
havior at small acceleration even though diB^ ^ 0. We 
choose the functional form of 11 (really valid only for the 
acceleration due to a point mass) as in @ 



/' 



however, now 



((V$) 2 -(^') 2 ) 1 



1-p 



(18) 



(19) 



In principle, p is an arbitrary function of K, but we take 
it to be a constant. Note that 



(V$) 2 > {^f 



(20) 



so the presence of an inhomogeneous vector field estab- 
lishes a minimum acceleration. Moreover, the closer p is 
to 1, the greater the influence of the vector field. The 
usual MOND-type solution is recovered if 



($B J ') 2 =P(V3) S 



(21) 



We will take this condition as a starting point. Ulti- 
mately, one must demonstrate that the vector-field dy- 
namics leads to this relation, at least approximately. 
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METHOD 

Outside a single point mass the gravity in the deep 
MOND regime is 



V Gjvmao « 
V $ = r 



(22) 



The solution of the vector field for equation (f2"Tj) is then 



B = ± pGjfmao 



(23) 



Equation (|23)l is not well defined at r = 0; however, since 
we are not going to discuss the acceleration near the cen- 
ter, we ignore this divergence. 

Of more concern is that the vector field A a should 
approach (1,0,0,0), in other words B — ► 0, at spatial 
infinity. It seems that this requires that either the grav- 
ity deviates from the MOND limit (equation [22]) at very 
large r, or p(r) — ► at large r. At any rate, we shall 
assume that MOND applies on cluster scales, and that p 
remains constant on such scales, so that any deviations 
from these two approximations occurs on scales much 
larger than the size of a cluster. 

As has been shown in Q, perturbations of the vector 
field grow with the cosmological scale factor. Structure 
is formed as the vector field perturbations attract matter 
particles and thereby cause the gravitational potential 
to grow. The size and the mass of a galaxy or a clus- 
ter depends on the size and amplitude of the vector field 
perturbation. The fate of the seed vector perturbations, 
however, is unstudied. If these spatial terms of the vector 
field do not decay after the cluster structure is formed, 
then the space surrounding the cluster should still have 
a vector field that is different than the value that would 
be inferred by solving the vector field equation with the 
matter as a source. We assume that the vector field per- 
turbation extends over a region that is larger than the 
size of our simulation, about 8Mpc. 

For two point particles with mass mi and m.2, the grav- 
itational acceleration in the MOND regime superposes 
as: 



(v*) 2 ~ (V*i) 2 + (v$ 2 ) 2 = 



(24) 

Note that this is very different than the usual Newtonian 
superposition of V*£- To satisfy this relation we assume 
that the vector fields adding according to: 

(d t BJ) 2 = (diB{) 2 + {diB{f, 

where B\ and Bi are the vector fields associated with 
(but not necessarily sourced by) m\ and m; . 

There are several numerical methods 17, li|, 2(| that 
can be used to calculate the gravitational acceleration 
in the environments of a mass distribution according to 



MOND. To simplify this problem, we choose an axially 
symmetric mass distribution. The method in 17[ is more 
suitable and easier for this case. The code solves the 
following equations: 



Here 



y -U = AttGnP 
V x vU = 0. 



U = v * 



(25) 
(26) 



(27) 
(28) 



serve as definitions of U and v respecitively. Here, /i is 
defined by equations (fl"8|) and (Tj"9|) . On a cylindrical 



grid, defined by {(rj, Zj)} (where i indexes the radial co- 
ordinate, r, and j indexes the axial coordinate z, and 
the symmetry allows us to ignore the azimuthal coordi- 
nate), these differential equations can be replaced by the 
following difference equations: 



(U, 



-i u i+i,j - r i U i,j) =- 4 ' Y '" 



-2it5zj(ri. 



(29) 



= 0. 



(30) 



Here U[j is the r component of U at 7"j, Zj\ similarly £7f 



is the z component. The r)k are discrete representations 
of fxdl, eg. rji = (i/j_i5rj_i + Vij8ri)j2. For details of 
the calculational method, please refer to [l7T |. 

Since our goal is to investigate whether the weak lens- 
ing maps of the bullet cluster can be reproduced in the 
MOND-like limit of GEA, we must calculate those maps. 
The lensing convergence is 



DdsDd 
c 2 D s 



V 2 <S>dl 



(31) 



Here, Dd, Dd s and D s are the angular distances between 
the observer and the lens, the lens and the source, and the 
observer and the source, respectively, i is the distance 
along the light ray from the observer. Once we have 
solved equations (|29p and (130)) , we can straightforwardly 
calculate k. 



MODEL AND RESULT 

Our code is able to simulate only axially symmetric 
mass distributions. The mass distribution of the bullet 
cluster, fortunately, is very close to axially symmetric. 
The angle between the direction of the motion of the 
bullet cluster and the the plane of the sky is about 8°. 
We, therefore, assume that the bullet cluster is axially 
symmetric about an axis that is parallel to the plane of 
the sky. 



5 





FIG. 2: The covergence map changing with different p: — 50kpc for the vector fields, p = 0.00, p — 0.50, p = 0.9, and 
p = 0.99 are corresponding to left up, right up, left bottom, and right bottom. The interval between each contour is 0.15. 
Left up figure shows the result of MOND. It does not show a gravitational lensing center at — 350kpc and 350fcpc. As the p 
increases, concentracted contours appears. 



For the baryonic mass distribution of the cluster, we 
take a King profile [28[ : 



p(r) 



Po 



i + (£) 5 



(32) 



We choose a core radius, r c = 150kpc. The central den- 
sity, po, and the positions of mass centers are chosen 
acording to table Q] as extracted from 21] . Figure Q] We 
choose l.hMpc to be the cutoff of the mass distribution. 
Dd £ ) Dd = 540-Mpc, which is appropriate for the lensed 
galaxies in the background of the bullet cluster, where the 
source of the gravitational lensing is at Zi — 0.296, and 
the a mean redshift of the lensed galaxies is Zb g = 0.85 

In addition to the baryonic clusters, we add vector field 
configurations at position A and D. The only coupling 
constants in the model are the usual gravitational cou- 
pling Gn, and the very small acceleration aQ. We there- 
fore assume that the interaction between the vector field 



"halos" associated with each cluster are sufficiently small 
that the two vector-halos can pass through each other 
unimpeded, much as dark matter halos are assumed to. 
The mass term in equation (|22p for the vector field at 
A is chosen to be the mass of component A plus that 
of component B. Similarly, the mass term in equation 
(f2"2"]l for the vector field at D is chosen to be the mass 
of component C plus that of component D. We assume 
the mass also follows equation ([3"2"]l . But the r^s of the 
vector fields are different from those of the baryonic mass 
distributions. 

Figure [5] shows the convergence map changing as p 
is varied. Here, we choose = bOkpc for the vector 
fields. For p = 0, the solution is the same as MOND 
There are two lensing centers that appear in the map, 
but they are located approximately halfway between A 
and B, and halfway between C and D. Thus, to repro- 
duce the lensing map of the bullet cluster would require 
vector-concentrations that are very far outside the bary- 



6 





FIG. 3: The convergence map changing with different for the vector fields: p = 0.99 for the model, — 200kpc, — lOOfcpc, 
th = 25kpc, and ru = Wkpc are corresponding to left up, right up, left bottom, and right bottom. The interval between each 
contour is 0.15. As the decreases, concentracted contours appears near — 350fcpc and 350£;pc. 



TABLE I: The mass and position of baryon matter 



component 


type 


position 
(kpc) 


Plasma mass 
within 100/spc 


A 


lensing 
center 


-350 


6.0 xlO 12 M 


B 


X-ray 
center 


-150 


6.8 xl0 12 M© 


C 


X-ray 
center 


200 


3.3 xl0 12 M Q 


D 


lensing 
center 


350 


5.9 xlO 12 M 



onic mass distributions. 

However, as p is increased, the centers of the conver- 
gence contours move away from the baryonic centers, and 
k is much larger than in MOND. This means that for p 
near 1, the vector fields play a role very much like dark 



matter, and are likely to augment the lensing in very 
much the same way as do dark matter halos. 

Figure [3] shows the convergence map changing as Th 
for the vector fields is varied. To emphasize the effect of 
the vector field, we choose p = 0.99. Moving from panel 
a (top left) through panel d (bottom right), changes 
from 200fcpc to lOkpc. For large r^, k is less than in 
MOND. This is because (from equation [T9f the effective 
ao is 

a e ff = a Q ^l -p- (33) 

The MOND acceleration regime applies only when a < 
a e ff. If the scale height of the vector field (ru) is very 
large, then a < a e ff only at very large distances; at 
smaller distances, the gravity becomes Newtonian, k is 
then much less than it would have been under MOND. As 
r/j is reduced, a gravitational center appears and moves 
to the centers of the vector fields and the n also becomes 
bigger. Clearly, Figure [5] and [3] suggest that the space 
components of vector fields can replace dark matter. 
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CONCLUSION 

In this work, we analyzed the effects of gravitational 
lcnsing caused by the vector field in the Generalized 
Einstein- Aether theory. The results show that the vec- 
tor field could cause the dissociation of mass center and 
gravity center observed in the bullet cluster. This would 
not be the vector field sourced by the matter, but rather 
the vector field that seeded the growth of the matter 
perturbation. We have not followed the evolution of this 
primordial seed, but rather established the approximate 
properties that it would have to explain the observations. 

The problem with the vector field is that a particular 
configuration of the vector field distribution is needed. 
This is the same problem as the need for a particular mass 
distribution of the dark matter. Future work will address 
comparing the full weak lensing data instead of the con- 
vergence map only. We note that, since TeVeS also in- 
cludes the same type of vector field, it is reasonable to 
believe that the gravitational acceleration of TeVeS could 
be dissociated from mass centers by the vector fields. 

The convergence map alone is not enough to distin- 
guish this type of vector field from dark matter. We 
need a more detailed analysis to distinguish the two mod- 
els, eg. by looking at the lensing shear we may yet find 
that modified gravity is distinguishable from dark mat- 
ter. Finally we note that this independent concentration 
of vector is very much like a dark matter halo - but a 
halo of vector field rather than of particles. Thus MOND 
started as an effort to build a theory in which what you 
see (baryons) is what you get (i.e. sources all gravity). 
One way or another, even in modified gravity theories, 
that seems not to be true. 

It is a pleasure to thank Constantinos Skordis, Irit 
Maor, Pedro Ferreira and Thomas Zlosnik for very useful 
discussions. We thank the authors of 2l| for permission 
to reproduce figure [TJ This work has been supported by 
grants from the US DOE, NASA and HEPCOS group at 
SUNY. 
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